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Abstract 



Presence of higher derivative terms in the Horava model of gravity can generate an 
instability in the Minkowski ground state. This in turn leads to a space dependent vacuum 
metric with a length scale determined by the higher derivative coupling coefficient. The 
translation invariance is spontaneously broken in the process. The phenomenon is interpreted 
as a form of Landau liquid-solid phase translation. The (metric) condensate acts as a source 
that modifies the Newtonian potential below the length scale but keeps it unchanged for 
sufficiently large distance. 



Introduction and Formalism : Horava Gravity (HG) [T] has created a lot of recent 
interest but surprisingly an intriguing consequence of the all important higher derivative 
extension terms in HG has gone unnoticed: possibility of an instability in the conventional 
flat metric ground state leading to a richer structure of the vacuum. This is a form of 
the celebrated Landau liquid-solid phase transition [2] of the spacetime (only space in the 
present context) itself. In the present work we show that this is indeed feasible and provide 

a preliminary study of some of its immediate effects. 
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An important lesson that we have learnt from Horava's [T] novel proposal of a UV com- 
plete model of Gravity is that Lorentz invariance is not that sacred in a theory. (See [3] 
for a lucid exposition of the perspective regarding Lorentz invariance violation.) It is ap- 
parently admissible even for a fundamental theory of Gravity, such as the Horava model, 
to be Lorentz non-invariant at ultra high energy or short distance, at ~ trans-Planckian 
regime, so long as Lorentz invariance and Einstein General Relativity (GR) is recovered at 
observable scales of low energy (or long distance). This has brought about a paradigm shift 
in the way of thinking of High Energy Physicists but is quite conventional for Condensed 
Matter Physicists where maintenance of Lorentz invariance is not an issue at all. Indeed, 
it was quite well known that higher derivative covariant terms, (generically functions of R, 
the Ricci scalar), in Einstein gravity improves the UV behavior but brings in insurmount- 
able ghost problems [3]. To overcome this in Horava Gravity (HG) only higher order spatial 
derivative terms are kept. These explicitly Lorentz breaking terms do not introduce ghosts 
and at the same time can cure the UV divergence problems. However the jury is still out 
on the question of whether HG is the true theory of gravity since the loss of Lorentz and 
diffeomorphism symmetries induces a spurious degree of freedom beside the spin 2 graviton 
[5]. Several ways to improve the Horava model have also been suggested [6]. 

On the other hand, the higher derivative nature of the Horava model opens up a com- 
pletely new line of thought that is of interest to us: instability in the flat metric ground 
state. This can lead to a Spontaneous Symmetry Breaking (SSB) as regards to translation 
symmetry only. This phenomenon is a liquid-solid type of phase transition in the space(time) 
itself a la Landau J^. A short distance length scale is generated in the process of transition 
from a homogeneous (~ liquid like) phase ground state to an inhomogeneous (~ solid like) 
condensate phase ground state in spacetime. Hence our work strengthens the claim of HG 
as a viable model for Quantum Gravity because the physics below this scale is affected while 
GR is recovered for distances sufficiently above this scale. The more ambitious project is 
to construct a true spacetime crystal where translation as well as rotation symmetries are 
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lost. The obvious motivation is to build a suitable spacetime for Quantum Gravity which is 
expected to have a short distance (~ Planck) scale □. Fortunately we know explicitly how 
to proceed in Condensed Matter Physics where the physics is naturally non-relativistic. 
Since we are in the Horava framework of gravity where Lorentz invariance is explicitly broken, 
formalisms exploited in (non-relativistic) Condensed Matter Physics should be applicable in 
the present scenario as well. We mention an earlier work by us p2] where a length scale 
was introduced in the metric in a higher derivative Lorentz covariant Gravity theory via 
SSB. However, the new non-trivial ground state broke the Lorentz invariance explicitly and 
furthermore the parent model has the well known ghost problem. In this sense HG serves 
perfectly as the laboratory to test this new idea since HG is free of ghosts and is not a 
Lorentz invariant theory to begin with (recall that the higher derivatives are only in the 
spatial sector). 

In a series of pioneering works Alexander and Mctague [11] were able to construct a 
solid crystalline lattice (where translation and rotation symmetries are lost), from a liquid 
phase (where the symmetries are intact), through SSB. This analysis explained the fact that 
BCC lattice structure is favored for metals undergoing liquid-solid transition, essentially in 
a model independent way. In order to achieve this the condensate has to have additional 
structure, i. e. its VEV can not be a constant that we generally encounter in Particle Physics. 
Rather, the crucial feature of the condensate is that the Fourier transform of non-zero VEV 
of condensate minimizing the free energy must have support at a non-vanishing momentum. 
Evidently, as is explained below, this requirement translates, in the coordinate space, to the 
parent lagrangian undergoing SSB, as having higher (at least fourth) order derivative terms 



3 The attempts so far made are some what ad-hoc: (i) Non-Commutative (NC) geometry framework [7] 
that exploits Seiberg Witten map [5] and incorporates NC corrections on GR. 

(ii) A cutoff length scale in gravity was introduced [9] by mapping the slice of spacetime into the phase space 
of quantum mechanics of fermions. The length scale and noncommutativity parameter in gravity map to 
the Planck's constant in quantum mechanics. 

(iii) A granular spacetime is considered in |10) . 
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that are quadratic in the field. The field in question is identified with the difference between 
the (inhomogeneous) density for solid and (constant) density of liquid, that acts as the order 
parameter. The free energy is expanded around the higher symmetry liquid phase and it 
should be mentioned that to form a proper crystal lattice one needs terms, at least of third 
and fourth order terms in the order parameter, in the expansion. However, the latter does 
not concern us in the present work as we will restrict ourselves to higher derivative quadratic 
terms in the HG action. The ideas developed in [TT] were applied in High Energy Physics in 
the context of string compactification by [13J. 

SSB of Higher Derivative Scalar Theory: As a warm up exercise we consider a toy 
model consisting of a scalar field theory with higher spatial derivative terms, 

s = Jd'x ^0(-«9 o 2 + d*)<p - ^mn - Vifi) 

In fact one can think of (CQ) as obtained from 

r n a 2 
s = Jd'x [-<p(-d 2 + d 2 )<p - ^(-d 2 + d 2 ) 2 <p - v{4>) 

by simply leaving out the time derivatives in a 2 term. This is in the same spirit as obtain- 
ing the Horava model with higher order spatial Ricci scalar terms from a covariant higher 
derivative gravity model containing higher orders of the full Ricci scalar. Furthermore, as 
we will see later, after some simplifying (but not unconventional) restriction on the metric 
fluctuation, the Horava model will be structurally identical to this toy model. In the present 
preliminary analysis we drop the potential V(4>) as does not play any role since we are con- 
cerned only with ground state solution obtained by minimizing the kinetic term. Minimizing 
V((f>) restricts the solution, obtained by minimization of the kinetic energy, still further. 
Explicit structure of V(4>), such as with cubic and higher order terms, becomes crucial in 
building the lattice in momentum space of (see [H] for more details). 

We are working in flat Minkowski spacetime with the metric rj^ = diag(—l, 1, 1, 1). The 
energy is given by 

H = Jd?x [\^ 2 + i(^0) 2 + (2) 
4 



As explained above we have dropped V(4>) from further consideration. 

Let us now minimize the energy. It is natural to consider the ground state to be static. 
Using the Fourier decomposition 

<P(x) = J^j <l>(k)e tk - x dk, 

the ground state energy can be written in the momentum space as 

1 - - a 2 -* 

h = --em^-k) + -{eym^-k). (3) 

Clearly H[k 2 = 0] = but H[k 2 = 1/a 2 ] = —4^2 showing that the minimum of kinetic 
energy fl3]) occurs at k 2 = \. Here we used the notation \k\ 2 = k 2 . 

From Figure 1, one can clearly see that the higher derivative term shifts the position of 
the energy minimum from k 2 = to k 2 = 1/a 2 . Subsequently translational invariance is 
broken in the higher derivative theory [TT1 [T3] . 



Mk 2 ) 




Figure 1: Plot of A(k 2 ) = k 2 vs k and A(k 2 ) = —k 2 (l — \k 2 ) vs k where the energies 
are respectively given by H(k 2 ) = <f)(k)<f)(—k)A(k 2 ) and H(k 2 ) = <p(k)<j)(—k)A(k 2 ). The 
specific form of the function A(k 2 ) is determined by the underlying microscopic theory [13]. 
Translational symmetry is present in the first case (dashed line) as the minima of A(k 2 ) 
occurs at k 2 = 0. Translational symmetry is broken in the second case (thick line) as the 
minima of A(k 2 ) occurs at k 2 ^ 0. 
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The condensate value of < <f)(x) >, (that minimizes the energy), has a nontrivial coordi- 
nate dependence which can be determined as follows: 



where C turns out to be a dimensionless number. In the above expression Q, notice that 
although the condensate <f>(\ x |) has become explicitly x-dependent, indicating that the 
translation invariance is lost, the space rotational symmetry is preserved. This is because in 
the Fourier integral the restriction was to integrate on the spherical surface k 2 = 1/a 2 that 
does not affect rotational symmetry. For large a energy H in ([3]) is minimum for k 2 = and 
the condensate < 0(| x |) > in Q becomes | x (-independent constant. 

It is customary to consider fluctuations above the condensate by shifting the field. In the 
present instance this exercise will means <p(x) — > 4>{x) — 0(| x |) leading to a translation 
non-invariant (but rotation invariant) theory. There is an underlying periodic nature of the 
translation symmetry broken phase. The coupling coefficient a in the higher derivative term 
in the action is responsible for the inhomogeneity and is present in the resulting lagrangian. 
We have not provided an explicit expression for that. Hence we have succeeded in introducing 
a periodic structure in the deformed translation invariance through SSB and the periodicity 
scale can be directly linked with the coupling coefficient of the higher derivative term. We 
now repeat the same exercise on the Horava model of Gravity in weak field approximation. 

SSB of Higher Derivative Horava Gravity: As we have explained above, breaking 
of translation invariance via SSB requires higher derivative terms in the action and this 
motivates us to consider higher derivative gravity theories. We take the higher derivative 
terms as perturbations on the linearized Horava gravity action (which reproduces Einstein's 






C sin(\ x | I a) 



(4) 



(2vr) 2 
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GR theory in the low-energy (IR) limit). As discussed before, Lorentz symmetry is already 
explicitly broken in Horava theory. 

Our long-term goal is to generate a spacetime lattice through SSB in the way described 
in the previous section. However, in the present article our target is more modest: as a 
first step towards this objective, we intend to introduce a short distance scale (~ Planck 
length) in the Horava action that will break the translation invariance, without affecting the 
rotational symmetry. Indeed, the extension of general relativity to NC spacetimes [7] also 
aims to achieve that but, compared to this ad- hoc procedure, our approach is much more 
basic and intuitive. We introduce the deformation directly in the metric as a fluctuation 
and the condensate plays the role of a (tensorial) order parameter in the continuum-discrete 
spacetime phase transition. As we have extensively discussed before, this requires the in- 
troduction of higher derivative terms in the action, R^R^ and (g^Rij) 2 . Apart from the 
conventional (g^-Rj^-term, the above are needed to ensure that the kinetic term is minimized 
for a non-zero momentum. 

We start with the action S 

S = J dtC = J dtd 3 Xy/gN(G ijkl K i:i K kl + AR + BRijR ij + CR 2 ) 

= J dtdfxy/gNiKijKV - XK 2 + AR + BR i:j R ij + CR 2 ). (5) 

Here is the spatial metric, A, B, C are three dimension-full parameters of the theory, 
N is the Lapse function, R is the spatial Ricci scalar and is the extrinsic curvature 
defined as 

Kij = 5ft®>9v ~ ViNj - VyJVi). (6) 

with Ni(x,t) is the Shift vector in ADM formalism [H] and the generalized De Witt metric 
Qijki j g defined as 

G ijkl =~(g ik g jl +g il g jk )-\g ij g kl . (7) 
For A = 1, A = 1 and B = C = Horava gravity reduces to Einstein's gravity. 
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Throughout we will work in the weak field approximation with 

g%j = Sij + hij , N = 1 + n , Ni = m. (8) 

Under these perturbations (jHJ), the expressions for the extrinsic curvature and the Ricci 
scalar turn out to be 



Kij = ^(dohij - diiij - djTii) , K = 8 13 K i:j = -(d h - 2d,, 



J, 



R {j = ^(d k dih jk + d k djh ik - d 2 hij - didjh) , R = didjh ij - d 2 h. (9) 
Using the above expressions Q and the relation 

y/gR = h tJ + ±S*R) (10) 

in the action Q we obtain the lagrangian density C of second order in h: 

C = h[dohijdoh* j - \(d h) 2 - 4(d n i )(6> J /i ij - X&h) + (4A - 2)n i (d i d j n j ) - 2md 2 n 11 



+jh ij (d 2 h lJ -2d k d i h jk +2d i &> h-5 ij d 2 h)+An(d i d J h i] -d 2 h)+C(d i d j h i] -d 2 h)(d k d l h M -d 2 h) 
+-{d k d i h jk - d 2 h ij + d k d 3 h ik - didjh) {d^ 1 h jl - d 2 h ij + d { &h il - &&h). (11) 

In the above expression Q and throughout the rest of our paper, we used the following 
notation: 

h = b i3 h i3 , d 2 = did 1 = 5 ij did r 
From this lagrangian Q we obtain the conjugate momenta as 

' s mrr° ■ ' fs e^rr-^ hV -^ (12) 
^^s(i^=5^-w- (13) 

Taking the trace of the relation ( II 3ft . we can write d^h^ in terms of tt^ as 

d h ij = 2 (n* + . (14) 



Using the relations (|T2|) . ( IT3|) . f|T4|) in Q we get the Hamiltonian density, 

U=p i (d Q n i )+T? i (d Q h ij )-C 

= - — — tt 2 - -(2A - l)m(&&nj) + -n;9V 

J 3A - 1 2 V ' v j; 2 

-\ j (d 2 ti j -2d k d i h jk +2d i d j h-5 ij d 2 h)-An(d i d j h ij - 

-j{d k dih jk - 9 2 /^- + d k djh ik - didjfyidi&h? 1 - d 2 h i] + di&W 1 - tf&h). (15) 
Using the momentum constraints and the gauge n.j = 0, the Hamiltonian Q now reduces to 

-j(d k dih jk - d 2 h i3 + 9*9^ - did ^{30 b? 1 - d 2 h i] + 9*9^ - 9*9^). (16) 

In the static limit 7Ty = and with /iy = t^/i for the sake of convenience, the Hamiltonian 
Q can be written in the Fourier space as % s tatic{p)i 

n static ( P ) = -jh(p)h(- P ) (i - ™ P 2 ) . (it) 

The minimum of this ( I17|) occurs at 

p 2 - 



a 2 MB' 

In a similar way as in Q we have the explicit form of f(r) as 

r 

where r = |x| stands for the radial coordinate. 

Now, to consider the SSB effects, let us consider fluctuations h^ v above the condensate, 

hij = Sijf(\S\) + hj{x) 

hoo = h 00 = -2<p , h 0i = 0, (19) 

where /(| x |) is a purely spatial function. The /(| x |) term, which is exactly of the same 
form as in Q, comes from the energy minimization arguments (previously described in the 



scalar theory case) since we do not want the presence of a constant vector in the theory as 
it would break the rotational invariance, and h^ u is treated as the order parameter. Q 
In terms of the metric Q, the lagrangian Q becomes 

£ = ^[dji^dji 11 - X(d h) 2 ] + jh l3 {d 2 ~ti j - 2d k d r h ]k + 2d l d j h - 5 ij d 2 h) 

+An(d i d j h ij - d 2 h) + Cididfti - d 2 h){d k d l h kl - d 2 h) 

+j(d k d i h jk - d 2 h io - + d k d J h i k - didfh^did'y 1 - d 2 h lj + dt&h* - &&h) 

- kd 2 f) ~ tChPididjtff - 5 l3 (d 2 ) 2 f) - —fCHfrdjPf - 5 t3 {d 2 ) 2 f). (20) 

Varying this lagrangian Q with respect to the fields n and we get the following equations 
of motion 

didft* - d 2 h - 2d 2 f = 0, (21) 



h&dohij ~ Xd d h) + ^(d% 3 - 2d k dih jk + didfh + 5 i:j d k &h k i ~ S tJ d 2 h) 
+A{d l d j n - Sij&n) + 2C(d l d j d k di~h kl - d 2 d l dfh - 8 ij d 2 d k d l h kl + S^fh) 

+j(2d i d j d k d l h kl - d 2 d k d-h 3k - d 2 d%h ik + (d 2 )%, - S ij ^d k d l h kl + d^fh) 

A 3 R 

+-(8,0, f - 5 l3 d 2 f) - ACid^d 2 } - 5 tJ (d 2 ) 2 f) - —{d&cPf - 5 lJ (d 2 ) 2 f) = 0. (22) 

Note that the /-terms in Q and in (12"T|) are subsequently acting as the components of a 
conserved source. 

Following standard literature [15J, we decompose h^ u as: 

h 00 = -2n = -20 , h i = , hij = 2s i3 - 2i/;S ij , (23) 



4 It is now clear that to make a full fledged spacetime "crystal" one requires third and fourth terms in 
huvi ( as explained in [111113] for liquid-solid transition), in the action which will considerably complicate the 
model. We are not considering this additional terms in the present work. 
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where Sij is traceless. Using the above decomposition (|23|) and considering the static limit 
(hence dropping all time derivatives), the equations of motion ( )2T| ^ turn out to be: 

d 2 ^ = \d 2 f, (24) 



Am - w - *) + A d % d + * m - 5 ^)f + *c#m - mv 

-ACd 2 (d % 8 3 - 5 tJ d 2 )f + B(d 2 ) 2 s l3 - 3B5 lJ (d 2 ) 2 i; - ^-d 2 (d l d J - 5 tJ d 2 )f = 0. (25) 



The solution for tp is obtained from ff24|) as 

GM l gm(|) 

where the first term in the right hand side of (1261) is the usual term as obtained in linearized 
Einstein Gravity and the second term is the modification due to the condensate. We now 
take the trace of both sides of the equation Q to obtain the following relation between and 

d 2 i<P-i> + t + 3 ^d 2 f] = 

<f> = 4 ~ - ~ ^rd 2 f (27) 
2 4A v ' 

Using the explicit expression of /(r) as given in (fl8|) . the expression for turns out to be 

GM 3B sinji) 

(p[r) = h — 7—z — . (28) 

r AAa 2 r 

The modified Newtonian potential derived above in (J2SJ) is shown in Figure 2. One can see 
that the fluctuations die out for sufficiently large distance r » a but can have non-trivial 
effects for r ~ a. 

Conclusion: To conclude, we have shown that higher derivative terms in the Ho- 
rava model of gravity can lead to an instability and a spontaneous breaking of translation 
symmetry of the space takes place. Hence ground state gets modified from the conventional 
Minkowski vacuum metric to a space dependent form. This is an example of Landau type of 
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Figure 2: P/ot o/ the usual Newtonian potential (thin line in the plot) and the modified 
Newtonian potential as described in previous paragraph (dashed line in the plot) vs r with 
a = 2. 

liquid-solid phase transition occurring in the space itself. A length scale also appears in the 
resulting metric in a natural way. The Newtonian potential is recovered for distances large 
compared to this scale. Consequences of this novel form of spatial metric is worth pursuing. 
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